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Abstract 

The dynamic dipole polarizabilities for the Li atom and the Be"*" ion in the 2 and 2 states 
are calculated using the variational method with a Hylleraas basis. The present polarizabilities 
represent the definitive values in the non-relativistic limit. Corrections due to relativistic effects are 
also estimated. Analytic representations of the polarizabilities for frequency ranges encompassing 
the n = 3 excitations are presented. The recommended polarizabilities for ''Li and ^Be^ were 
164.11 ± 0.03 al and 24.489 ± 0.004 ag. 

PACS numbers: 31.15.ac, 31.15.ap, 34.20.Cf 
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I. INTRODUCTION 



The advent of cold atom physics has lead to increased importance being given to the pre- 
cise determination atomic polarizabilities and related quantities. One very important source 
of systematic error in the new generation of atomic frequency standards is the blackbody 
radiation (BBR) shift [l]-[3[. The differential Stark shifts caused by the ambient electro- 
magnetic field leads to a temperature dependent shift in the transition frequency of the two 
states involved in the clock transition. The dynamic polarizability is also useful in the deter- 
mination of the magic wavelength in optical lattices 0^ • Another area where polarization 
phenomena is important is in the determination of global potential surfaces for diatomic 
molecules jsj. 

When consideration is given to all the atoms and ions commonly used in cold atom 
physics, the Li atom and Be"*" ion have the advantage that they have only three electrons. 
This makes them accessible to calculations using correlated basis sets with the consequence 
that many properties of these systems can be computed to a high degree of precision. The 
results of these first principle calculations can serve as atomic based standards for quantities 
that are not amenable to precision measurement. For example, cold-atom interferometry 
has been used to measure the ground state polarizabilities of the Li and Na atoms [9|, 
lOj. However the polarizability ratio, ad{X) / adiLi) can be measured to a higher degree 
of precision than the individual polarizabilities [ll]. So measurements of this ratio, in 
conjunction with a high precision ab-initio calculation could lead to a new level of accuracy 
in polarizability measurements for the atomic species most commonly used in cold-atom 
physics. 

Calculations and measurements of Stark shifts are particularly important in atomic clock 
research since the BBR shift is predominantly determined by the Stark shift of the two levels 
involved in the clock transition. The best experimental measurements of the Stark shift have 
been carried out for the alkali atoms and accuracies better than 0.1% have been reported [l2|, 



13| . Experimental work at this level of accuracy relies on a very precise determination of the 



electric field strength in the interaction region |13l-ll5j . High precision Hylleraas calculations 



of the type presented here provide an invaluable test of the experimental reliability since 
they provide an independent means for the calibration of electric fields (l6| . 

The dynamic Stark shift in oscillating electromagnetic fields is also of interest. The so- 
called magic wavelength, i.e. the precise wavelength at which the Stark shifts for upper 
and lower levels of the clock transition are the same, is an important parameter for optical 
lattices. The present calculation is used to estimate the magic wavelength for the Li 2 "^S 
— 7- 2 transition. The present calculations of the AC Stark shift potentially provides an 
atomic based standard of electromagnetic (EM) field intensity for finite frequency radiation. 

There have been many calculations of the static polarizabilities of the ground and excited 
states of the Li atom and the Be"*" ion The most precise calculations on Li and Be^ 

are the Hylleraas calculations by Tang and collaborators [l7|, |23| . The Hylleraas calculations 
were non-relativistic and also included finite mass effects for Li. Large scale calculations 
using fully correlated Hylleraas basis sets can attain a degree of precision not possible for 
calculations based on orbital basis sets 19|, |2J, |25| . There have been many calculations of the 



dynamic polarizability for Li [isl, 126-32 , but fewer for Be^ His. The present calculation is 
by far the most precise calculation of the dynamic polarizability that is based upon a solution 
of the non-relativistic Schrodinger equation. One particularly noteworthy treatment is the 
relativistic single-double all-order many body perturbation theory calculation (MBPT-SD) 
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TABLE I: Comparisons of the binding energies (in a.u.) of Li and Be"*" in their low-lying states. 
The experimental valence binding energies are taken from the National Institute of Standards 
database [sij]- The J- weighted average is used for states with L > 1. The ground-state en- 
ergies for the °°Li+ and °°Be'^+ ions are -7.2799134126693059 and -13.6555662384235867 a.u. 
respectively [3j]. The ground-state energies for ^Li"*" and ^Be^"*" are —7.2793215198156744 and 
—13.6547092682827917 a.u. respectively [3j]. Underlining is used to indicate digits that have not 
converged with respect to basis set enlargement. 



State 




Theory 




Experiment 




<^Li 




^Li 


6'7Li 


2^5 


-0.19814691124 




-0.19813041084 


-0.198142 


2 2p 


-0.13024311963 




-0.13023623876 


-0.130236 


325 


-0.07418381350 




-0.07417777025 


-0.074182 


3^P 


-0.05723769823 




-0.05723424577 


-0.057236 


3^D 


-0.05561012974 




-0.05560578543 


-0.055606 


425 


-0.03752870957 




-0.03752445073 


-0.038615 




-0 03139073613 




-0.03138814390 


-0.031975 


4 2d 


-0.03127588444 




-0.03127343938 


-0.031274 


42^ 


-0.03125355531 




-0.03125111202 


-0.031243 




°°Be+ 




9Be+ 


9Be+ 




-0.66919693847 




-0.66915422599 


-0.669247 


2^P 


-0.52376705352 




-0.52375065365 


-0.523769 


325 


-0.26720549176 




-0.26718867334 


-0.267233 


3 2p 


-0.22956788615 




-0.22955822005 


-0.229582 


3^D 


-0.22248781972 




-0.22247429085 


-0.222478 


425 


-0.13629487843 




-0.13628082370 


-0.143152 


42p 


-0.12222924451 




-0.12221999823 


-0.128134 


42^ 


-0.12512688879 




-0.12511926908 


-0.125124 


42F 


-0.12501546711 




-0.12500785769 


-0.125008 



by Safronova et al. [32[. This calculation is fully relativistic and treats correlation effects 
to a high level of accuracy, although it does not achieve the same level of precision as the 
present Hylleraas calculation. 

The present work computes the dynamic dipole polarizabilities of the Li atom and the 
Be+ ion in the 2 25, and 2 2p levels using a large variational calculation with a Hylleraas basis 
set. This methodology allows for the determination of the computational uncertainty related 
to the convergence of the basis set. Analytic representations of the dynamic polarizabilities 
are made so they can subsequently be computed at any frequency. Finally, the difference 
between the calculated and experimental binding energies is used to estimate the size of the 
relativistic correction to the polarizability. The final polarizabilities should be regarded as 
the recommended polarizabilities for comparison with experiment. All quantities given in 
this work are reported in atomic units except where indicated otherwise. 
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II. THE STRUCTURE CALCULATIONS 



A. Hamiltonian and Hylleraas coordinates 

The Li atom and Be"*" ion are four-body Coulomb systems. After separating the center 
of mass coordinates, the nonrelativistic Hamiltonian can be written in the form SS"! 



--Ei;;v?-i-gv..v,-tf 

^ 1 



i>j>l "-^ 

where r^j = |rj — rj\ is the distance between electrons i and j, fi = mome/(mo + is the 
reduced mass between the electron and the nucleus, and Z is the nuclear charge. In our 
calculation the wave functions are expanded in terms of the explicitly correlated basis set in 
Hylleraas coordinates: 

A^fr- r r- \ — r^'^ ^jl2 ^j23 ^331 „-ari- I3r2-tr3 

(P\'-li '-2, l3) — '1 '2 '3 '12 '23 '31 ^ 



X 



where 3^(^^^)^^2 is the vector-coupled product of spherical harmonics to form an eigenstate 
of total angular momentum L and component 

^{^SKi2A(^l'^2,r3) = ^ (^l"^l4?«2|^12"^12) 

all rrii 

X {iurriuisms \ LML)Ye^rm (ri) Y^^m^ {r2)Ye^rn3 (h) , 

(3) 

and ^(l, 2, 3) is the three-electron spin 1/2 wave function. The variational wave function is 
a linear combination of anti-symmetrized basis functions 0. With some truncations to avoid 
potential numerical linear dependence, all terms in Eq. (|2]) are included such that 

3l + 32 + J3 + jl2 + j23 + 731 < ^ , (4) 

where Vt is an integer. The computational details in evaluating the necessary matrix elements 
of the Hamiltonian may be found in j25|. The nonlinear parameters a, /3, and 7 in Eq. ([2]) 
are optimized using Newton's method. 

The convergence for the energies and other expectation values is studied by increasing Vl 
progressively. The basis sets are essentially the same as two earlier Hylleraas calculations of 



the static polarizabilities [17|, |23[. The maximum f2 used in the present calculations is 12. 
The uncertainty in the final value of any quantity is usually estimated to be equal to the 
size of the extrapolation from the largest explicit calculation. 

Fig. [U is a schematic diagram showing the nonrelativistic energy levels of the most im- 
portant states of the Li atom. The energy level diagram for the low lying states of Be"*" is 
similar. 

The energies of the ground states for °°Li and ^Li were —7.47806032391(5) and 
—7.47745193065(5) a.u. respectively. The respective energies for the "^Be"*" and ^Be"*" ground 
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FIG. 1: Low lying energy levels of the Li atom. The energy level diagram for Be"*" is similar. 



states were -14.3247631769(3) and -14.3238634942(3) a.u.. Table [T] gives the binding en- 
ergies of the Li atom and Be"*" ion systems with respect to the two-electron Li+ and Be^"*" 
cores. The Hylleraas basis was optimized to compute the 2 and 2 state polarizabilities, 
so some of the n = 4 state energies have significant deviations from the experimental n = 4 
state energies. The states with significant energy differences can be regarded as pseudo- 
states. The uncertainties listed in Table [T] represent the uncertainties in the energy with 
respect to an infinite basis calculation. The actual computational uncertainty is very small 
and there is no computational error in any of the calculated digits listed in Table [H 

With one exception, all the finite mass binding energies are less tightly bound than 
experiment. The differences from experiment are most likely due to relativistic effects. The 
exception where experiment is less tightly bound than the finite mass calculation is the 4 '^F 
state of Li. This exception was not investigated since the properties of this state do not 
enter into any of the polarizability calculations. 

B. Polarizability definitions 

The dynamic polarizability provides a measure of the reaction of an atom to an external 
electromagnetic field. The dynamic polarizability at real frequencies can be expressed in 
terms of a sum over all intermediate states, including the continuum. The dynamic dipole 
polarizability is expressed in terms of the dynamic scalar and tensor dipole polarizabilities, 
ai{uj) and aj{uj), which can be expressed in terms of the reduced matrix elements of the 
dipole transition operator: 




(5) 



L. 



'a 




(6) 



L. 



'a 



where 



ai{La,uj) 



9(2L + 1) 



n 



AEoJ{noL\\Ti\\nLa) 




(7) 
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with Ti = X]i=o 9'*-^*^io(R-i) being the dipole transition operator, and 

W{L,La) = 




30(2L+ 1)L(2L - 1^ 



(2L + 3)(L + 1) 

In the above, \nQL) is the initial state with principal quantum number no, angular momen- 
tum quantum number L, and energy Eq. The nth intermediate eigenfunction |nLa), with 
principal quantum number n and angular momentum quantum number L^, has an energy 
En- The transition energy is Aii^on = En — Eq. The qt are the charges of the respective 



particles and are defined in Ref. [35!]. In particular, for the case of L = 



ai{uj) = ai(P,w), (9) 

afM = 0; (10) 

for L = 1, 

ai{uj) = ai{S,u) + ai{P,u) + ai{D,u) , (11) 

a^{u) = -ai{S,uj) + ^ai{P,uj) - ^a,{D,u). (12) 

In Eqs. (fTTj) and f|T2l) . ai{P, u) is the contribution from the even-parity configuration [pp')P. 
The scalar and tensor polarizabilities can be easily related to the polarizabilities of the 
magnetic sub- levels, ai^M^^^), 

aifi{u) = ai{u) — 2aj{u) 
ai,±i{co) = ai{u) + ai{u) . (13) 

III. THE DYNAMIC POLARIZABILITY FOR THE ~LI ATOM AND THE ~BE+ 
ION 

A. Ground state dynamic polarizabilities 

Fig. [2] shows the dynamic dipole polarizability of the lithium ground state as a function of 
photon energy. The chief errors in the dynamic polarizability are related to the convergence 
of the n^P excited state energies. The largest calculation used a basis with dimensions 
{Ns,Np) = (6412,5761). The difference between the ai^u) and polarizability computed 
with a {Ns,Np) = (4172,3543) basis would be barely discernible in Fig. [2l The convergence 
of ai{u) is best at photon energies far from the discrete excitation energies of the n^P 
excitations. The polarizability is very susceptible to small changes in the physical energies 
at photon energies close to the n^P excitation energies. 

The uncertainties in the dynamic dipole polarizabilities of the Li ground state as well as 
the polarizabilities themselves are listed in Table [III All of the values listed are accurate to 
about ±1 in the fifth digit for u < 0.11388 a.u.. Some of the alternate calculations of the 
ai{u) polarizabilities 18, 261-28^32 [are listed in Table [ITl Dynamic polarizabilities from 



some less accurate calculations "Bill have not been tabulated 
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0.00 0.05 0.10 0.15 

(0 (a.u.) 



FIG. 2: Dynamic dipole polarizability, ai(cj), of the Li atom in the ground state. The singularities 
in the polarizabiUty at the 2^5 ^ n^P frequencies are marked. 
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FIG. 3: The dynamic dipole polarizability, a\{bS) for the ground state of the Be"'' ion. The singu- 
larities in the polarizability at the 2 ^5" — ?> n frequencies are marked. 



One feature of Table [IT] is the excellent agreement with the MBPT-SD calculation of 
Safronova et al. 33|- The MBPT-SD calculation and the present Hylleraas calculation are 
in perfect agreement when the MBPT-SD theoretical uncertainty is taken into consideration. 
While the MBPT-SD calculation is fully relativistic, its treatment of electron correlation is 
less exact than the present calculation. The MBPT-SD calculation also gives no consider- 
ation of finite mass effects. Relativistic effects would tend to decrease ai{u) at low u, and 
the MBPT-SD calculation gives slightly smaller ai{u) at low u. 

The older CI-Hylleraas calculation of values of Pipin and Bishop j26| compares excel- 
lently with the present more modern calculation. All digits in 0:1(0;) from the CI-Hylleraas 
calculation are in perfect agreement with the present Hylleraas calculation. The model po- 
tential polarizabilities of Cohen and Themelis [18| are also very close to the present dynamic 
polarizability. The Cohen-Themelis potential was constructed using a Rydberg-Klein-Rees 
(RKR) inversion method. The Time-Dependent-Gauge-Invariant (TDGI) polarizabilities of 
Merawa et al [23, 28| are only accurate to 0.5% or larger. The moderate accuracy of TDGI 
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TABLE II: The dynamic dipole polarizabilities, ai{uj) (in a.u.), for the Li ground state. The results 
of the fourth column include relativistic corrections. The numbers in brackets for the second and 
third columns are the uncertainties in the last digits arising from incomplete convergence of the 
basis set. The uncertainties in the recommended (Rec.) values reflect additional uncertainties 



related to the relativistic correction. 







Hylleraas 




MBPT-SD [32] TDGI 


CI-Hyllcraas 


Model 




°°Li 




Rec. ^Li 


[27, 2^ 


L26] 


Potential [18] 


0.00000 


164.112(1) 


164.161(1) 


164.11(3) 


163.6 


164.1 


164.14 


0.00500 


164.996(1) 


165.045(1) 


165.00(3) 


164.5 


165.0 


165.03 


0.01000 


167.707(1) 


167.758(1) 


167.71(3) 


167.2 


167.7 


167.74 


0.02000 


179.517(1) 


179.574(1) 


179.52(3) 


178.9 


179.5 


179.55 


0.02931 


201.242(2) 


201.313(1) 


201.24(3) 


201.0(7) 






0.03000 


203.438(1) 


203.512(1) 


203.44(3) 


202.6 


203.4 


203.47 


0.03420 


219.221(1) 


219.307(1) 


219.22(4) 


219.0(8) 






0.04000 


250.265(1) 


250.376(1) 


250.26(4) 


248.8 


250.3 


250.29 


0.04624 


304.278(1) 


304.441(1) 


304.26(5) 


304.0(8) 






0.05000 


356.077(1) 


356.300(1) 


356.05(6) 


355.2 


356.1 


356.60 


0.05699 


550.259(1) 


550.790(1) 


550.18(9) 


549.7(1.1) 






0.06000 


741.165(2) 


742.126(1) 


741.00(12) 


729.2 




740.73 


0.06507 


1984.577(1) 


1991.488(1) 


1983.11(31) 


1983(3) 






0.07000 


-2581.603(2) 


-2569.994(2) 


-2584.54(40) 


-2895.3 






0.07592 


-645.478(2) 


-644.749(1) 


-645.70(10) 


-645.9(1.3) 






0.08000 


-415.067(1) 


-414.763(1) 


-415.17(7) 


-427.1 






0.09000 


-211.518(2) 


-211.439(1) 


-211.55(3) 


-216.5 






0.09110 


-199.941(1) 


-199.868(1) 


-199.97(3) 


-200.1(0.8) 






0.10000 


-135.872(2) 


-135.838(1) 


-135.89(3) 


-0.819 






0.11388 


-86.266(1) 


-86.249(1) 


-86.27(2) 


-86.4(0.8) 






0.15183 


-38.210(9) 


-38.204(9) 


-38.22(1) 


-38.4(1.1) 






0.16000 


-31.08(5) 


-31.06(5) 


-31.08(6) 









calculations has also been noted in calculations of the static polarizabilities [17|. 

The static polarizabilities for Be"*" in the infinite mass approximation have been presented 
recently [23[. The present calculation represents an extension of this earlier calculation 
since finite mass effects are now included. The dynamic dipole polarizabilities listed in 
Table IIIII includes transition frequencies that extend well into the ultraviolet region. The 
most accurate of the few alternate calculations should be the CI-Hylleraas calculation of 
Muszynska et al [30[. However, it gives an ai{uj) that is about 1% smaller than the present 
polarizability. Space limitations precluded tabulation of the TDGI polarizability [28|. The 
TDGI polarizability was of only moderate accuracy with errors of about 2% for u < 0.6 
a.u.. The uncertainty in the present polarizability is about 10~^ a.u. for photon energies 
lower than 0.40 a.u., but has increased to 10^^ a.u. at a; = 0.50 a.u.. 

The dynamic polarizabilities for the Li and Be"^ ground states are depicted in Figures 
|2] and |3l There are obvious similarities in shapes of the two ai [u) curves but with the Li 
polarizability being about 5-10 times larger in magnitude at comparable values of uj/u}2s^2p- 
One difference between the two curves is that Be"*" has zeroes in ai{u) at a discernible 
frequency difference before the 3 and 4 excitations while the ai {u) negative to positive 
crossovers for Li occur much closer to the transition frequencies. 
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TABLE III: Dynamic dipole polarizabilities, ai{uj) (in a.u.), for the ground state of the Be+ 
ion. The results of the fourth column incorporate relativistic effects. The numbers in brackets 
are the uncertainties in the last digits arising from incomplete convergence of the basis set. The 
recommended (Rec.) polarizabilities in the fourth column reflect uncertainties other than purely 
computational. 



UJ 




Hylleraas 




CI-Hylleraas 


(a.u.) 


°°Be+ 


9Be+ 


Rec. ^Be+ 


[301 


0.00 


24.4966(1) 


24.5064(1) 


24.489(4) 


24.3 


0.01 


24.6088(1) 


24.6187(1) 


24.601(4) 


24.4 


0.02 


24.9518(1) 


24.9620(1) 


24.943(4) 


24.7 


0.04 


26.4291(1) 


26 4404fl) 


26 419(4) 


26.2 


0.06 


29.3390(1) 


29.3528(1) 


29.325(5) 


29.1 


0.08 


34.7358(1) 


34.7550(1) 


34.715(6) 


34.3 


0.10 


45.6509(1) 


45.6836(1) 


45.609(7) 


44.9 


0.12 


74.7857(1) 


74.8724(1) 


74.656(12) 




0.15 


-367.8708(2) 


-365.8030(3) 


-371.860(60) 




0.18 


-43.2038(1) 


-43.1745(1) 


-43.273(7) 




0.20 


-25.3195(1) 


-25.3090(1) 


-25.348(4) 




0.30 


-5.7967(1) 


-5.7951(1) 


-5.801(2) 




0.40 


-0.2912(1) 


-0.2873(1) 


-0.2961(2) 




0.50 


-2.149(7) 


-2.161(7) 


-2.164(8) 





B. Excited state dynamic polarizabilities 

The scalar and tensor dipole polarizabilities for the excited 2 state of the Li atom 
are listed in Table IIVI As far as we know the present calculations are the only dynamic 
polarizabilities presented for this state. The structure of the dynamic polarizability is com- 
plicated since both downward and upward transitions leads to singularities. This is seen 
most clearly in Fig. |4] which plots the polarizabilities for photon energies up to 0.10 a.u.. 
The tensor polarizability is generally small except in the vicinity of the 2 ^5, 3 "^S and 3 
transitions. The tensor polarizability can become large when a single transition tends to 
dominate Eq. f lT^ . The scalar and tensor polarizabilities tend to be opposite in sign. The 
main contribution to the polarizabilities comes from transitions to the S and D states. The 
coefficients in the sum- rules, Eqs. (fTT!) and (fT2!) . for these terms are opposite in sign. 

The dynamic polarizabilities for the Be"*" 2 state are also tabulated in Table [IV] and 
depicted in Figure [5] for the photon frequencies below 0.40 a.u.. There are three resonances 
in this frequency range. The scalar and tensor dynamic polarizabilities are similar in shape 
but with the opposite sign. As far as we know, there has been no previous calculation of 
the 2 ^P state dynamic polarizability. 
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FIG. 4: The dynamic polarizabilities, ai(w) and af{Lo) (in a.u.) of the Li 2^P state for photon 
frequencies below 0.10 a.u.. The scalar polarizability is given by the solid line while the tensor 
polarizability is given by the chain curve. 




FIG. 5: The dynamic polarizabilities, ai{u}) and aj{uj) (in a.u.), of the Be"^ 2^P state for photon 
frequencies below 0.40 a.u.. The scalar polarizability is given by the solid line while the tensor 
polarizability is given by the chain curve. 

C. The static 2 25' ^ 2 2p Stark shift 

The static Stark shift for the 2 ^S* — 2 energy interval in an electric field of strength 
F is written as 

A-E2s-2p,M = ~2^^ ^"^^s ~ a2p,M) 

- (72s - 72p,m) + . . . , (14) 

where 7 is the hyper-polarizability. The Stark shift depends on the magnetic quantum 
number M of the 2 state. The relative size of Aa and A7 determines the extent to which 
the Stark shift is influenced by the hyper-polarizability at high field strengths. The relative 
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TABLE IV: The dynamic dipole polarizabilities of the 2^P state of Li and Be^. Both the scalar 
and tensor polarizabihties are tabulated. The numbers in brackets are the uncertainties in the last 
digits arising from incomplete convergence of the basis set. Values without uncertainties have no 
numerical uncertainties in any of the quoted digits. The recommended (Rec.) polarizabilities in the 
sixth and seventh columns have estimated corrections from relativistic effects. The recommended 
polarizabilities reflect uncertainties other than purely computational. 





°°Li 


^Li 




Rec. 


^Li 


(a.u.) 




aT 


ai 


aT 


ai 


aT 


0.00 


126.9458(3) 


1.6214(3) 


126.9472(5) 


1.6351(2) 


126.970(4) 


1.612(4) 


0.01 


129.2491(5) 


1.4035(2) 


129.2501(5) 


1.4178(2) 


129.273(4) 


1.393(4) 


0.02 


136.8371(5) 


0.5302(3) 


136.8372(5) 


0.5463(5) 


136.864(4) 


0.518(4) 


0.03 


152.469(1) 


-2.091(1) 


152.468(2) 


-2.070(1) 


152.503(4) 


-2.106(4) 


0.04 


185.542(5) 


-11.722(5) 


185.535(5) 


-11.691(5) 


185.593(6) 


-11.747(6) 


0.05 


301.24(8) 


-82.33(9) 


301.23(9) 


-82.27(8) 


301.33(10) 


-82.38(10) 


0.06 


-119.1(5) 


446.7(5) 


-119.3(5) 


446.9(5) 


-119.0(6) 


446.6(6) 


0.07 


1804.5(1) 


-904.2(2) 


1801.2(1) 


-900.34(5) 


1806.2(2) 


-905.2(2) 


0.08 


-593.1(2) 


-43.5(4) 


-592.7(3) 


-43.6(5) 


-592.6(5) 


-43.4(5) 




oo 


Be+ 


9Be+ 


Rec. 


9Be+ 




ai 


aT 


Oil 


aT 


ai 


aT 


0.00 


2.02476(1) 


5.856012(1) 


2.02319(1) 


5.858938(1) 


2.0285(10) 


5.8528(10) 


0.01 


1.99755(1) 


5.890887(1) 


1.99595(1) 


5.893842(1) 


2.0013(10) 


5.8876(10) 


0.02 


1.91389(1) 


5.997630(1) 


1.91221(1) 


6.000672(1) 


1.9178(12) 


5.9942(12) 


0.05 


1.23144(1) 


6.845876(1) 


1.22905(1) 


6.849643(1) 


1.2363(13) 


6.8415(13) 


0.10 


-3.88178(1) 


12.61790 


-3.89080(1) 


12.62842 


-3.8666(15) 


12.6033(13) 


0.14 


-94.71454 


104.48595 


-95.2476 


105.02073 


-93.7873(16) 


103.5592(16) 


0.20 


26.1505(1) 


-12.9481(2) 


26.1499(1) 


-12.9451(1) 


26.161(17) 


-12.958(12) 


0.25 


48.59(5) 


-26.10(1) 


48.61(5) 


-26.11(1) 


48.62(6) 


-26.17(1) 


0.28 


52.09(1) 


-3.43(1) 


52.10(1) 


-3.44(1) 


52.07(1) 


-3.45(1) 


0.32 


-49.87(1) 


4.413(1) 


-49.85(1) 


4.412(1) 


-49.89(1) 


4.40(1) 



importance of Aa and A7 is given by the ratio 

X = = Elp , (15) 

12{a2s - a2p,M) 12Aq; 

Using the static polarizability and static hyper-polarizability for the Li atom results in 
Aa = 37.1 and A7 = 9.99 x 10^ giving X = 0.0001 at F = 6.67 x 10"^ a.u. (344 kV/cm) 
and X = 0.001 at F = 2.11 x 10"^ a.u. (1087 kV/cm). These estimates of the critical field 
strength where the quadratic Stark shift is valid depend slightly on the magnetic quantum 
number and exact values can be determined by using M-dependent polarizabilities. Stark 
shifts of higher order than the hyper-polarizability can be comfortably ignored at the 0.01% 
level provided the field strength is less than 1100 kV/cm. The static Stark shift for Be"*" is 
not interesting since it is difficult to measure as a Be"*" ion immersed in a finite electric field 
is accelerated away from the finite field region. 
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FIG. 6: The polarizability difference between the 2 '^S and 2 states of Li. Polarizabihty differences 
are shown for M = and M = 1. 
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FIG. 7: The polarizabihty difference between the 2 and 2 states of Be"^ 
differences are shown for M = and M = 1. 



Polarizabihty 



D. The dynamic 2^5 ^ 2 Stark shift 

The Li Stark shifts, a{2s) — a{2pM), are plotted as a function of frequency in Figure El 
It is seen that there are magic wavelengths for M = just below the 2 — )■ 3 ^S* threshold 
and between the 2 ^5 — )■ 2^P and 2^P — )■ 3^D thresholds. The actual energies for which 
the polarizabihty difference is zero are given in Table |Vl The Stark shifts get very large for 
frequencies between 0.058 and 0.070 a.u.. 

The Be+ Stark shifts, a{2s) — a{2pM), are plotted as a function of frequency in Figure [71 
The Stark shifts are much smaller in magnitude than the Li atom shifts. One difference 
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TABLE V: The photon energies for which there is no Stark shift for the — ?> 2^P transition. 
UnderHned digits indicate uncertain digits arising from lack of basis set convergence. Digits in 
brackets indicate possible uncertainties associated with relativistic corrections in the recommended 
(Rec.) values. 



System 



M = 



M = 1 



^Li 



^Li 



Rec. ^Li 



^Be+ 



Rec. ^Be+ 



0.0463176806 
0.0810217955 
0.0936643305 
0.0463356878 
0.0810244789 
0.0936618991 
0.046297(4) 
0.081014(2) 
0.0936613(2) 
0.2629202678 
0.3703715027 
0.3907521463 
0.2629173603 
0.3702799522 
0.3904553568 
0.2628956(7) 
0.370274(1) 
0.3904546(2) 



0.0847639571 



0.0847660870 



0.084756(2) 



0.3784570004 



0.3784518437 



0.378443(2) 



from Li is that the Be"*" shift has no zero for energies below the 2 ^S* — )■ 2 threshold. The 
first zero in the Stark shift (excepting those related to a singularity) is at 0.263 a.u.. 

E. Analytic representation 

The utility of the present calculations can be increased by constructing a closed form 
expression for the dynamic polarizability. This is done by retaining the first 3 terms in Eq. ([7]) 
explicitly and then expanding the energy denominator in the remainder. The expressions 
explicitly include oscillator strengths up to the n = 4 principal quantum numbers. The 
closed form expression is 

+ uj^S{-Q) + ...+ uj^^S{-lQ) + C{u) (16) 

where 

^(-m) = y-%^^, (17) 
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TABLE VI: The parameters for the calculation of the 2 '^S state frequency-dependent polarizabilities 
of Li and Be"*". The numbers in the square brackets denote powers of 10. The recommended (Rec.) 
results of the fourth and the seventh columns incorporate relativistic corrections. 



Parameter 


°°Li 




Ree. ■^Li 


°°Be+ 


9Be+ 


Rec. 9Be+ 


/2s->2p 


0.746956855381 


0.746961871867 


0.747011776131 


0.498067422721 


0.498083382699 


0.498227010322 


A-E2s2p 


0.067903791567 


0.067894172078 


0.06790605 


0.145429884364 


0.145403572344 


0.14547806 


/2s->3p 


0.004731019443 


0.004737600312 


0.004728028090 


0.083243986131 


0.0832889414647 


0.083209271939 


A-E2s3p 


0.140909212964 


0.140896165068 


0.14090640 


0.439629051730 


0.439596005937 


0.43966521 


J2s—i-4:p 








u.u^uo 1 ^uooyui 






AE2s4p 


0.166756175058 


0.166742266938 




0.546967693380 


0.546934227760 




S(-2) 


1.69771 


1.69863 




0.379627 


0.379781 




S(-4) 


21.8714 


21.8888 




0.513938 


0.514207 




S(-6) 


428.809 


429.237 




0.979476 


0.980088 




S(-8) 


9.69364(3] 


9.70502(3] 




2.05322 


2.05469 




S{~W) 


2.37008(5] 


2.37325(5] 




4.52145 


4.52508 




S(-12) 


6.08098(6] 


6.09006(6] 




10.2473 


10.2564 




S(-14) 


1.61032(8] 


1.61296(8] 




23.6368 


23.6596 




^(-le) 


4.35760(9] 


4.36537(9] 




55.1280 


55.1852 






27.0605 


27.0643 




2.33229 


2.33246 





dco) = . (18) 



1 — r/iw^ 

Here f2s^np are the dipole oscillator strengths for the 2 ^S* — )■ n transitions with transition 
energies AE2snp- The S{—n) are the Cauchy moments of the remainder of the oscillator 
strength distribution and are independent of u. The C{uj) is an approximate term to rep- 
resent the summation from the term S'(— 18) to ^(oo). The ratio, rji = S{~n — 2)/S{—n), 
is assumed to be constant and its value is set at S'(— 16)/5'(— 14). Numerical values of the 
various constants in Eq. f[T^ can be found in Table IVI[ Inclusion of the remainder term has 
greatly increased the precision of the analytic fit to the exact dynamic polarizability. 

The analytic representation for the Li state is accurate to 0.01 a.u. for u ^ 0.1612 
a.u. and to an accuracy of 0.1 a.u. for u ^ 0.1728 a.u.. The dynamic polarizability for the 
Be"*" 2^5* state maintains its accuracy over a larger u range. It is accurate to 0.001 a.u. for 
u ^ 0.543 a.u., to 0.01 a.u. for oj ^ 0.58605 a.u. and 0.1 a.u. for u ^ 0.6 a.u.. 

The presence of zeroes in the dynamic polarizability near the singularities means that the 
relative error in the analytic representation can get very large in a frequency range very close 
to the zeroes. Neglecting these localized regions with anomalously high relative uncertainties, 
the relative difference between the analytic representation and actual dynamic polarizability 
for the Li 2^5 state was less than 0.001% for to ^ 0.1399 a.u., 0.01% for co ^ 0.1551 a.u., 
and 0.1% for u ^ 0.1651 a.u.. The relative difference for the Be"*" 2^5* state obtained by the 
variational Hylleraas method was less than 0.001% for u ^ 0.4737 a.u., 0.01% for u ^ 0.5067 
a.u. and 0.1% for u ^ 0.52575 a.u.. The inclusion of the remainder term, C{u) improved 
the accuracy of the analytic representation by one or two order of magnitude within the 
frequency range listed above. 

The dynamic dipole polarizabilities of the 2 states of Li and Be"*" have both scalar. 
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ai{u), and tensor, af{u), parts. The scalar part can be written 



/ \ \ ^ f2p^ns , \ ^ f2p- 



where 



+ ^(-2) + io^S{-4) + uj^S{-6) + ... 

+ io^^S{-lQ) + C{uj), (19) 



n=5 



The 2p — )■ n'P excitation involves a core excitation and the intermediate state is an 
unnatural parity ^P^ state. The tensor part is 

4 4 

X/ \ /2p— s-ns 1 f2p-^nd 

+ w^^S^(-16) + C^H, (21) 



where 



cTf \ _ \ ^ f2p^ns 1 \ ^ f2p ^n'P 

b [ m) - 2^ _-—_- + - 2^ __ 



„=5 (AE2p„,)- 2 ^ (AE2p„'p)- 
J_ \ ^ f2p->-nd 

10 ^ (AB,^,)» • ^ ' 



(23, 

where f2p^mLi means the oscillator strength from 2p state to mLi state transition. 5''^(— 2), 
S'"^(— 4), S'^(— 6) ■ ■ ■ are the coefficients corresponding to a;", w^, ■ ■ ■ terms of the tensor 
part. The remainder term, C^iuj) is an approximate expression to take into account the 
S'^{—18) — )■ S'^(oo) summations. The factor is set to be T^f = S'(— 16)/S'(— 14). All 
parameters in the analytic representation are given in Table I VI II 

The ffist two terms of Eqs. f|T9|) and fl2T]) include five resonances, which make the major 
contribution to the polarizability with the second term involving excitations to D states 
being the most important. This is clearly seen in the Li ai{u) of 185.542(5) a.u. at w = 0.04 
a.u.. The contribution of the ffist summation of Eq. (1191) was —8.8717 a.u., while the second 
summation contributed 175.8241 a.u.. The value given by Eq. (1191) was 185.5378 a.u., which 
agrees with the exact value at the level of 0.0004%. 

The analytic representation for the scalar polarizability ai (u) of the Li 2 state is 
accurate to 0.01 a.u. for u ^ 0.0855 a.u. and to 0.1 a.u. for u ^ 0.0937 a.u.. The analytic 
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TABLE VII: The parameters defining the frequency-dependent polarizabilities of the 2 state of 
Li and Be"*". The numbers in the square brackets denote powers of 10. The recommended (Rec.) 
results of the fourth and the seventh columns incorporate relativistic effects. 



Parameter 


°°Li 




Rec. ''Li 


°°Be+ 


9Be+ 


Rec. 9Be+ 


f2p — 2 s 


-0.248985618454 


-0.248987290622 


-0.24900392538 


-0.166022474240 


-0.166027794233 


-0.166075670107 


A_E2p2s 


-0.067903791567 


-0.067894172078 


-0.067906050 


-0.145429884364 


-0.145403572344 


-0.145478060000 


/2p— )-3s 


0.110578835460 


0.110575403831 


0.110589306872 


0.064385095804 


0.064385347515 


0.064407168594 




0.056059306121 


0.056058468507 


0.056054150 


0.256561561765 


0.256561980308 


0.256536125000 


flp — ^4s 


0.014979087136 


0.014981586569 




0.012876689561 


0.012879003186 




A_E2p4s 


0.092714410055 


0.092711788029 




0.387472175091 


0.387469829948 




f2p—^3d 


0.638568044661 


0.638583007678 


0.638583083728 


0.631981700709 


0.632059480294 


0.632047210702 


0.074632989884 


0.074630453329 


0.074630150 


0.301279233806 


0.301276362793 


0.301291440000 


J ZifJ — r ^u. 


0.122746501135 


0.122756411337 




0.122711398708 


0.122734598764 




5(-2) 


0.098967235189 


0.098962799378 




0.398640164736 


0.398631384564 




16.8408 


16.8447 




1.07595 


1.07628 




5(-4) 


974.832 


975.135 




3.69533 


3.69658 




S(-6) 


6.40829 [4] 


6.41083(4] 




14.9422 


14.9479 




5(-8) 


4.48404[6] 


4.48620(6] 




64.3838 


64.4113 




'-'V ^"1 


O . J- W O 1 O J 






288 547 


288 684 




S( — 'l 2) 
'-'V ^^1 


2 43429 f1 01 






1 328 1 8 


1 328 87 




S(-14) 


1.85141 [12] 


1.85282(12] 




6232.21 


6235.79 




S(-16) 


1.42796[14] 


1.42918(14] 




29667.5 


29685.9 




m 


77.1282 


77.1354 




4.76034 


4.76058 




S^i-2) 


-2.73075 


-2.73118 




-0.156423 


-0.156454 




S^(-4) 


-154.670 


-154.702 




-0.544987 


-0.545110 




5^(-6) 


-9.88199[3] 


-9.88458(3] 




-2.10990 


-2.11044 




S^(-8) 


-6.71603(5] 


-6.71818(5] 




-8.71572 


-8.71824 




S^{-10) 


-4.73561(7] 


-4.73742(7] 




-37.5083 


-37.5206 




S^(-12) 


-3.42263(9] 


-3.42417(9] 




-166.164 


-166.226 




S^(-14) 


-2.51902(11] 


-2.52034(11] 




-752.681 


-753.003 




S^(-16) 


-1.88083(13] 


-1.88197(13] 




-3471.37 


-3473.06 




T 

m 


74.6651 


74.6713 




4.61201 


4.61228 





representation for the tensor polarizability, oc{{u), is accurate to 0.01 a.u. for u ^ 0.0926 
a.u. and to 0.10 a.u. for u ^ 0.1 a.u.. 

The relative error in the analytic representation of ai(cj) for the 2 state of the Li atom 
is less than 0.001% for uo ^ 0.082 a.u., 0.01% for u ^ 0.0871 a.u., and 0.1% for uj ^ 0.0906 
a.u.. The relative error of the analytic representation for ot{{u) is less than 0.001% for 
uj ^ 0.0815 a.u., 0.01% for uj ^ 0.0932 a.u., and 0.1% for u ^ 0.0995 a.u.. 

The dynamic polarizability of the Be+ 2 ^P state maintains its accuracy over a larger 
range of uj. It is accurate to 0.001 a.u. for uj ^ 0.3513 a.u., to 0.01 a.u. for uj ^ 0.3837 a.u. 
and 0.1 a.u. for u ^ 0.4111 a.u.. The absolute error for (w) is 0.001 a.u. for u ^ 0.3788 
a.u., 0.01 a.u. for a; ^ 0.4077 a.u., and 0.1 a.u. for uj ^ 0.4291 a.u.. 

The relative error between the analytic representation and Hylleraas values of ai{u) for 
the Be+ 2^P is less than 0.001% for u ^ 0.332 a.u., 0.01% for u ^ 0.3536 a.u. and 0.1% for 
io ^ 0.3708 a.u.. The relative error for oiJ{(jj) of the Be"*" 2^P state is less than 0.001% for 
UJ ^ 0.3265 a.u., 0.01% for uj ^ 0.3434 a.u. and 0.1% for uj ^ 0.3534 a.u.. 
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IV. FINITE MASS CORRECTIONS 



The effect of the finite mass was to decrease the Li atom and Be+ ion binding energies 
hsted in the Table HI Therefore, it is not surprisingly that the a; = polarizabilities of the 
Li and Be+ ground states are increased in Tables [III and lllli The overall changes of the 
a; = polarizabilities are 0.03% and 0.04% for Li and Be"*" respectively. The finite mass 
polarizabilities are larger than the infinite mass values at a; = a.u.. These differences can 
be taken as indicative of the overall change in the polarizabilities at finite frequencies below 
the first excitation threshold. The differences are naturally larger near thresholds. 

The finite mass effect for the Li 2^P state increased its polarizability by 0.001% (Table 
IIVP while decreasing the polarizability for the Be"*" 2^P state (Table [IV]) by 0.08%. This 
behavior for Be"*" is due to the 2^P — )■ downward transition. The increased negative 
contribution from this transition is enough to outweigh the increased positive contributions 
from transitions to more highly excited states. 

As a general rule, the magnitude of the polarizabilities for both upward and downward 
transitions increase for the finite mass calculations. The residual Cauchy moments S{—n) in 
Table |Vl] and Table IVIII all increase for the finite mass calculations since these are computed 
exclusively from upward transitions. 



V. OTHER EFFECTS AND UNCERTAINTIES 
A. Estimate of Relativistic effects 

The major omission from the present calculation is the inclusion of relativistic effects. 
The larger part of the energy difference between the present finite mass calculations and 
the experimental binding energies in Table |T] is due to the omission of relativistic effects. 
Relativistic effects will alter the polarizability calculation in two ways. First, the energy 
differences will be changed. Generally, the binding energies of all states can be expected 
to be slightly larger. Secondly, there will be some changes in the reduced matrix elements. 
The wave functions for the n "^S and n states can be expected to be slightly more compact 
since they are more tightly bound. 

Correcting for the relativistic energy is simply a matter of replacing the theoretical ener- 
gies in the sum rules by the experimental values. The spin-orbit weighted averages were used 
for states with L > 1. The corrections to the transition matrix elements are made by re- 
course to calculations using a semi-empirical model potential that supplements the potential 



field of a frozen Hartree-Fock (HF) core with a tunable polarization potential [20|, l23|, |36 



Polarizabilities for Li and Be"*" computed with this approach reproduce Hylleraas calculation 



at the 0.1% accuracy level [13, |20|, l23l . 

The method used to estimate the relativistic effect upon matrix elements relies on com- 
paring two very similar calculations. One calculation has its polarization potentials tuned 
to reproduce the finite mass energies of Table HI The other calculation is tuned to give 
the experimental energies. The matrix elements for the low lying transitions that dominate 
the dynamic polarizabilities are then compared. The differences between the "finite-mass" 
calculation and the "experimental" calculations are then determined. These changes in the 
matrix elements are then applied as corrections to the set of Hylleraas matrix elements. 
The only matrix elements that are changed are those involving transitions inside the 2 
and 3 level space. Transitions to these states dominate the 2 "^S and 2 ^P polarizabilities. 
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The actual change in the Li 2^5* —t- 2 matrix element was a reduction of 0.0054%. The 
reduction in the Be"*" 2^5" — ?■ 2^P matrix element was 0.011%. 

Using the new set of corrected matrix elements gives a ground state polarizability of 
164.114 a.u. (Table H]). This represents a reduction of the polarizability by 0.047 a.u.. A 
coupled cluster calculation of the lithium ground state estimated that relativistic effects 
reduced its polarizability by 0.06 a.u. 37|. The static polarizability of the 2^P state of ^Li, 
namely 126.947 a.u., was increased to 126.970 a.u. (Table HV|) . This gives a Stark shift of 
—37.144 a.u., which is in agreement with the experiment of Hunter et al jisj which gave 
—37.14(2) for the ^Li 2^5* — 2^Pi/2 Stark shift. Another calculation was made to check the 
2^A/2i2^P3/2 polarizability difference. The MBPT-SD calculation gave a difference of 0.015 
a.u. [22|. Doing two calculations tuned to give a 2 ^P spin-orbit splitting of 1.77 x 10~^ 
a.u.. (the energy splitting in the MBPT-SD calculation [22|) gave a polarizability difference 
of 0.0145 a.u.. A further test was made by examination of the line strengths of the Si^"*" 
3 ^S*- 3 ^P spin-orbit doublet. A MBPT-SD calculation gave a line strength ratio of 1.000524 
(once angular momentum factors were removed) [38[. Turning the core potential in a semi- 
empirical model based on the Schrodinger equation gave a value of 1.000618. The 
available evidence supports the conjecture that it is possible to use the energy differences 
between the Hylleraas and experimental energies to get an initial estimate of relativistic 
corrections for other properties such as the polarizability. The uncertainty in the correction 
would seem to be about 20%. To a certain extent the cancelations involved in adding the 
finite mass and relativistic corrections together leads to polarizabilities that are close to the 
infinite mass polarizabilities. 

The static polarizability of the Be^ ground state was reduced from 24.506 to 24.489 a.u.. 
This represents a reduction of 0.4%. However, the static scalar polarizability of the 2 ^P state 
increased from 2.0231 to 2.0285 a.u., an increase of 0.24%. The static tensor polarizability 
changed from 5.8589 to 5.8528 a.u.. The heavier mass and larger nuclear charge means 
relativistic effects are substantially larger than finite mass corrections. 

Dynamic polarizabilities and their analytic representations from the set of matrix ele- 
ments with the estimate of the relativistic effect are listed in Tables [III - IVHI as the rec- 
ommended values. The changes to analytic representation only involved changes in the 
oscillator strength and energy differences for a few states. 



TABLE VIII: Experimental C3 values from analysis of the Li2 spectrum and C3 values from the 
Hylleraas calculations. 



Source 



Value 



Li2 spectrum [39], Cg fixed from [19 
Li2 spectrum [8], Cq fixed from [19] 
Li2 spectrum [8], fixed from [20] 
Hylleraas °°Li 
Hylleraas ^Li 

Hylleraas ^Li: Recommended 



11.0022(24) 
11.00241(23) 
11.00240(23) 
11.000221 
11.001853 
11.0007 
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B. The 2^5* — ?• 2^P matrix element and uncertainties 

Recently Le Roy et al [8] analysed the ro-vibrational spectrum of the lithium dimer 
obtaining an estimate for the C3 parameter describing the long range C^/B? potential of 
the A-state that dissociates to the 2 "^S and 2 states. The C3 parameter can be related to 
the 2 ^5 — 2 multiplet strength. The determination of Le Roy et al represented an order 
of magnitude improvement in precision over any previous determination of C3. 

The current value of C3 = 11.0007 a.u. computed with relativistic corrections is about 
0.0155% smaller than the experimental value of Le Roy et al. The finite mass calculation 
with the relativistic correction is closer to experiment than the infinite mass C3, but there 
is a remaining discrepancy of 0.0017 a.u.. It is not likely that QED effects can explain the 
discrepancy as Pachucki et al found that these were 2.5 times smaller than relativistic effects 
in the polarizability of helium [40]. It must be recalled that the Le Roy et al experiment 
is reporting an order of magnitude improvement in experimental precision. Going to such 
extreme levels of precision means there might be small corrections that need to be applied to 
the analysis of the data that have not received consideration. For example, the value of C3 
will be different for states asymptotic to the 2 and 2 ^P3/2 levels. The analysis of Le Roy 
et al uses a common C3 value for both members of the spin-orbit doublet. Irrespective of 
this, it should be noted that experiment and theory are incompatible at precisions better 
than 0.01%. 

The difference between the present and Le Roy C3 is used to assign an error to the 
present polarizability calculation. Changes of 0.008% were made to the 2(3)^5 — 2(3) ^P 
matrix elements, the polarizabilities were recomputed, and the differences assigned as the 
uncertainty in the recommended values. This difference is actually larger than the estimated 
relativistic change in the matrix element. Therefore, the recommended static polarizability 
of the Li ground state is 164.11(3) a.u.. Uncertainties in the 2(3) — 3^D matrix elements 
are smaller (relativistic effects have a smaller impact on these matrix elements) and have not 
been included in the uncertainty analysis. The final value for the static scalar polarizability 
of the 2 2p state was 126.970(4) a.u., while the tensor polarizability was 1.612(4) a.u.. 

The same uncertainty analysis was applied to the Be"*" ion polarizabilities. The recom- 
mended value for the ground state is 24.489(4) a.u.. The static 2 ^P scalar polarizability was 
set as 2.0285(10) a.u. while the tensor polarizability was set to 5.8528(10) a.u.. 

Uncertainties in the recommended dynamic polarizabilities in Tables HTl IIIIl llVl andlVlwere 
computed by making corrections to the matrix elements, recomputing and then observing 
the change. These uncertainties are best interpreted as indicative as opposed to rigorous 
estimates. 

VI. SUMMARY 

Definitive non-relativistic values for the dynamic dipole polarizabilities of Li and Be+ in 
their low-lying 2 "^S and 2 ^P states have been established using the variational method with 
Hylleraas basis sets. Calculation for both finite and infinite nuclear mass systems have been 
performed. Analytic representations for the dynamic polarizabilities of the Li atom and the 
Be"*" ion have also been developed. These results can serve as a standard against which any 
other calculation can be judged. 

Subsidiary calculations have been used to estimate the impact of relativistic effects that 
are not explicitly included in the Hylleraas calculation. It is recommended that the value of 
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164.11(3) a.u. be adopted as the static polarizability of ^Li. The uncertainty of 0.03 a.u. is 
based on the difference between the present C3 and that of Le Roy et a/ [8|. This accuracy 
level is also supported by the ^Li 2 ^S" — 2 Stark shift of —37.14 a.u. which is in perfect 
agreement with the value of —37.14(2) given by the high precision experiment of Hunter et 



al [13|. The recommended static polarizability for Be^ is 24.489(4) a.u.. 

The dynamic polarizabilities that have been obtained can be used as an atom based stan- 
dard for electromagnetic field intensity. These polarizabilities can be regarded as an initial 
attempt to develop atom based standards for polarizability and Stark shift measurements. 
The primary virtue of the method with which the relativistic corrections were evaluated 
was simplicity of computation. For present purposes, the estimate of the relativistic correc- 
tions only has to be accurate to 10-20% for the recommended polarizabilities to be valid. 
The comparisons that have been done with fully relativistic calculations suggest that the 
estimates of the relativistic corrections are indeed accurate at this level. However, a more 
rigorous estimate using the Briet-Pauli Hamiltonian and perturbation theory would be de- 



sirable 40-42 
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